Universal mechanism of discontinuity of commensurate-incommensurate transitions in 
three-dimensional solids: Strain dependence of soliton self-energy 
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We show that there exists a universal mechanism of long-range soliton attraction in three- 
dimensional solids and, therefore, of discontinuity of any commensurate-incommensurate (C-IC) 
phase transition. This mechanism is due to the strain dependence of the soliton self-energy and 
specific features of the solid-state elasticity. The role of this mechanism is studied in detail for a 
class of C-IC transitions where the IC modulation is one-dimensional, the anisotropy in the order 
parameter space is small, and the symmetry of the systems allows the existence of the Lifshitz in- 
variant. Two other mechanisms of soliton attraction are operative here but the universal mechanism 
considered in this paper is found to be the most important one in some cases. Comparison with the 
most extensively studied C-IC transition in K2Se04 shows that the experimentally observed thermal 
anomalies can be understood as a result of the smearing of the theoretically predicted discontinuous 
transition. 
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PACS number: 64.70.Rh, 64.60.-i 

I. INTRODUCTION 

The earliest theories of commensurate-incommensurate 
(C-IC) transitions exposed a scenario of continuous 
transformation.^ According to this classical scenario, C- 
IC transitions occurred because the self-energy of some 
kind of soliton became negative so that these solitons ap- 
peared spontaneously. The interaction between soUtons 
was found to be completely repulsive. In consequence, 
the C-IC phase transition temperature coincided with 
the temperature of the sign change of the soliton self- 
energy, and the soliton density at this temperature was 
zero. 

The physical meaning of these solitons is different for 
different systems: they can be dislocations, domain walls, 
Abrikosov vortices, etc. The main result of the present 
paper is valid for all of them in three-dimensional solids, 
and it is deduced in detail considering these solitons as 
domain walls. This type of solitons appears, for in- 
stance, in IC dielectrics, IC magnetics, charge-density- 
wave systems and spin-Pierls compounds. Although the 
total number of these examples is difhcult to estimate, 
it numbers in the hundreds. What seems to be more 
clear is that the classical scenario has never been ob- 
served in C-IC transitions: the density of domain walls 
in IC phases is never small experimentally. Even in the 
so-called "domain-like regime" the distances between do- 
main walls are normally comparable with the domain wall 
widths and in the more frequently observed "sinusoidal 
regime" the domain walls strongly overlap and loose their 
individuality. 

Partially, this discrepancy between theoretical predic- 
tions and observations is purely experimental. It is 
because IC systems need a tremendously long time to 
achieve equilibrium. Its relaxation implies creation and 
arrangement of domain walls. The driving force of this 



arrangement is the interaction between the walls that, 
according to the above mentioned theories, decay expo- 
nentially with the interwall distance. At the same time 
the relaxation is expected to be strongly hampered due 
to pinning of the walls at defects and effects of discrete- 
ness of the crystal lattice. In other words, the theories 
deal with equilibrium states while in experiments one ob- 
serves, most probably, systems which are quite far from 
the equilibrium. 

It is also natural to suspect that some of the features 
of the classical scenario are results of implied approxima- 
tions and neglect of the theories. So one may ask what 
parts of this scenario will survive within a more complete 
theory, a question that arose long ago. There were many 
proposals to explain why the C-IC transitions should be 
discontinuous. There were too many, in fact, some of 
them proved to be unjustified. But still there are two sur- 
vivor effects, in some sense universal, which lead to the 
discontinuity: (i) the interaction between the order pa- 
rameter and the elastic strains via the strain dependence 
of the coefficient of the so-called Lifshitz invariant^''^ and 
(ii) the long-range attraction between the walls due to 
their thermal and quantum fluctuations.^ In this paper 
we show that a third effect must be taken into account. 
It is also a coupling between the order parameter and 
the elastic strains, but this coupling is described by the 
striction term. Any symmetry allows this coupling inde- 
pendently of the Lifshitz invariant allowance in the cor- 
responding Landau thermodynamic potential. We shall 
compare the effects of this striction-mediated mechanism 
and of the two mentioned above. The relative impor- 
tance of the three mechanisms is found to be different for 
different classes of C-IC transitions, with the striction- 
mediated mechanism the most important at least for 
some of them. 

We mention, first of all, that if the symmetries of the 
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two phases are not related the discontinuity of such a 
transitions is trivial. We suppose, therefore, that the 
structure of the IC phase can be presented as a space 
modulation of the structure of the commensurate phase, 
i.e., the two phases are describable in terms of the order 
parameter of the (implicit) normal-commensurate tran- 
sition. Then one has to distinguish between two cases, 
that were labelled as C-IC transitions of type I and type 
II in Ref.^. If the Lifshitz invariant is absent, IC-C tran- 
sitions (type II) are clearly discontinuous in the close- 
ness of the Lifshitz point. Michelson^ was the first who 
exposed this conclusion, which was confirmed later by 
many authors (see, e.g., Refs.^'^). However, under special 
conditions it has been predicted that type-II IC phases 
can obtain domain-like structures with continuous C-IC 
transitions. Recently such a behavior has been discussed 
for so-called spin-Pierls compounds. * In this paper we 
shall discuss in detail cases in which the Lifshitz invari- 
ant is present and the IC modulation appears along one 
crystallographic direction. The C-IC transitions in these 
cases (type I) are well established as continuous transi- 
tions if the above mentioned mechanisms are not taken 
into account. Systems that belong to the discussed 
cases number in the several tens of C-IC transitions of 
IC dielectrics^^ and there are at least several examples 
among magnetic^^ and charge density wave systems. 

The above described classical scenario has been mod- 
ified historically following the three Unes of the efforts 
clearly outlined in Ref.^: (i) overcome the constant- 
amplitude approximation, (ii) study the effects of fluc- 
tuations neglected within the Landau-like approach and 
(iii) analyze the role of coupling between the order pa- 
rameter and other variables. 

The development of studies along the first line was 
far from being smooth. Naive attempts to overcome the 
constant-amplitude approximation concluded that the 
transition should be discontinuous. However, subsequent 
studies showed that the performed procedure was incor- 
rect. The full story can be found in Ref.®. The most 
consistent and thorough efforts to overcome the constant- 
amplitude approximation were made by Golovko.^"'^^ His 
conclusion was, in short, that the C-IC transition remains 
continuous even well beyond the constant-amplitude ap- 
proximation in all cases he was able to treat exactly. 

Along the second line it was noted first that in the 
two-dimensional case fluctuations of domain walls en- 
hance the repulsive character of their interaction,^^ while 
these fluctuations are ineffective in the three-dimensional 
case.^® But it was found later that the combined effect of 
elastic and electric long-range fields and fluctuations of 
the walls results in a power-law attraction between the 
domain walls and, therefore, in the discontinuity of the 
C-IC transition.'* 

The studies along the third line was also not free of 
confusion. Bruce, Cowley, and Murray^ stated that the 
coupling of the order parameter to any "long wavelength 
coordinate" leaded to the discontinuity of C-IC transi- 
tions. The authors mentioned specifically the electri- 



cal polarization and the elastic strain, the latter via the 
strain dependence of the coefficient of the Lifshitz in- 
variant. Bak and Timonen-'^'^ criticized this result argu- 
ing that taking into account spatial changes of the ad- 
ditional coordinate one results in, simply, a renormal- 
ization of the coefficients of the thermodynamic poten- 
tial that does not change the character of the transition 
(Bruce, Cowley and Murray forced the coordinate to be 
uniform while allowed space variations of order param- 
eter). Paradoxically, Bak and Timonen referred to the 
elastic strain as this additional coordinate but their crit- 
icism is perfectly valid for any coordinate just excluding 
the elastic strain that need a special treatment.'''** The 
matter is in the long-range character of the elastic forces 
in solids according to the fact that to characterize elas- 
tic deformations of solids one uses different variables for 
the spatially homogeneous and spatially inhomogeneous 
parts of the deformations: six components of the strain 
tensor for the homogeneous part and three components 
of the displacement vector for the inhomogeneous one. 
Bak and Timonen overlooked the uniform part of the 
deformations while Bruce, Cowley and Murray excluded 
from consideration the inhomogeneous part. Subsequent 
works'''*^ overcome this controversy taking into account 
both types of deformations and showing that the results 
of Ref.^ remained qualitatively correct. 

If Bruce, Cowley and Murray had referred their the- 
ory to a solid with infinite shear modulus their results 
would be strictly correct. We shall use this possibility 
to illustrate the origin of the striction-mediated attrac- 
tion between solitons considered in the present paper. If 
the shear modulus of a solid is infinite its only possible 
deformation is an homogeneous dilatation e. Suppose 
that a finite density n of domain walls, or more generally 
solitons, is created in the solid. The soliton self-energy 
depends, naturally, on e. Taking the state without soli- 
tons as the non-deformed one, we present the change of 
the system energy per volume unit as 

F{u)c^{Eo + Eie)n+^e'', (1) 

were the terms in parentheses represent the soliton self- 
energy and K is the bulk modulus of the system. Mini- 
mizing F with respect to u one finds that the equilibrium 
deformation is Ceq = —nE\/K. Then, the change in en- 
ergy becomes into 

F{€eq)-Eon-^n^, (2) 
2K 

where the second term represents the attraction men- 
tioned before. ^° Evidently one can not forget that some 
soliton repulsion exists as well, a repulsion that provides 
a finite value of the equilibrium soliton density. Below we 
reproduce this soliton attraction, in cases in which these 
solitons represent domain walls, taking into account both 
homogeneous and inhomogeneous deformations as well as 
the elastic anisotropy of real crystals. 
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II. C-IC PHASE TRANSITION 



A. Elasticity-induced attractive interaction 



To describe the C-IC phase transition we shall use 

the Landau free energy expressed in terms of the or- 
der parameter of the "virtual" normal-commensurate (N- 
C) transition. We restrict ourselves to the case of two- 
components order parameter 77 = (p cos 772 = psinip) 
and an IC modulation along only one crystallographic 
direction. In a coordinate frame with the z axis along 
the IC modulation, the Landau thermodynamic poten- 
tial can be written as (see, e.g., Ref.^) 



F — Fi + F2, 



(3a) 



+5 



{^z^f + i^zpf] } dv, (3b) 



F2 = - f [rijUijp^ -I- gijUijp^Vzf + -^XijkiUijUki } dv. 



(3c) 



Here v is the volume of the system, a = ariT — Q) is 
the only temperature-dependent coefficient, 7 is related 
with the anisotropy in the order parameter space, and 
m > 3 is an integer that depends on the system in ques- 
tion. It is necessary to assume that the coefficients /3 
and S are positive from the stability conditions for the 
free energy. In Eq. (3c) Uij represents the strain tensor 
describing deformations of the normal phase (summation 
over double indices is implied). The term r-ijUijp^ takes 
into account the strain dependence of the coefficient q, 
i.e., the striction effect, and the term gijUij p^V zip the 
strain dependence of the coefficient of the Lifshitz invari- 
ant. 

Without these dependencies on the strain, the N-IC 
transition takes place at a temperature Ti > Q defined 
by the condition ao = axiT^ - 6) = (7^/(4:6). This N-IC 
transition is due to the presence of the Lifshitz invariant 
because if a were zero only a N-C transition should take 
place at T = ©. At the N-IC transition temperature Tj, 
the "wave number" of the IC structure is qq = a/ [25). 
Because of the anisotropy in the order parameter space, 
as the temperature decreases the IC structure becomes a 
domainlike one and finally, at some temperature Tloc^ the 
domain walls disappear in the IC-C transition (so-called 
lock-in transition). 

If we minimize Eq. (3) with respect to all elastic de- 
grees of freedom, the resulting thermodynamic potential 
contains nonlocal terms (see below), which makes it dif- 
ficult to determine the spatial distribution of the order 
parameter. But due to the specific features of an isotropic 
system with infinite shear modulus, to minimize Eq. (3) 
for such a system results quite simple having in mind al- 
ready known results. As we shall see, the results obtained 
for this "hypothetical" case can be easily generalized for 
real ones. 



1. Isotropic system with infinite shear modulus 

We consider first the case of an isotropic system with 
infinite shear modulus. Its only possible deformation is 
an homogeneous dilatation u. If we minimize Eq. (3) 
with respect to u we obtain 

(4) 

where Vij = rSij, gij = gdij, K represents the bulk mod- 
ulus of the system and (...) means volume average. As 
we have mentioned, further minimization of the Landau 
potential is not very easy due to the contribution of those 
nonlocal terms. 

But in this special case in which the shear modulus of 
the system is infinite there is another way of minimiza- 
tion. Note that the Landau potential Eq. (3) can be 
written as 

F =\ j {oi{'^)p^ + f^P'^ + 7P" cos{mip) - c7{u)p^VzV> 

'p'iVzvf + iVzpfWdv+^u^ (5) 
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where a{u) = a + ru, a{u) = a + gu, and u is a vari- 
ational parameter. Fixing for a moment this parameter 
one obtains the same equations of equilibrium for the 
order parameter as those, e.g., in Ref.^. Their solution 
gives, in particular, the thermodynamic potential close 
to the C-IC transition as a function of the domain wall 
density n. In the weak anisotropy case (see below) this 
can be written as^ 



K 

Fc{u) + —u^ + E{u)n- 



4:J{u)nex.p 



mp{u) 



2n 



(6) 



where the coefficients of this expression, in terms of 
the square of the order parameter amplitude p^{u) = 
—a{u)/{20) and the "wave number" of the IC structure 
at the N-IC transition go(w) = cr(u)/(2(5), are 



Jiu) 



^27 
5 



p'{u) = ^p"-^(u) 



E{u) 



Jiu) 



mqo{u) 

2p{u) - nqo{u)] 



2p{u) 



(7) 
(8) 

(9) 

(10) 



Let us analyze the physical meaning of Eq. (6) for a 
clamped crystal, i.e., a crystal with a fixed deformation 
u. The first two terms represent the free energy of the C 
phase (n = 0). The domain- wall self-energy is obtained 
from the third one. It is positive in the C phase, vanishes 
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at the transition point, and is negative in the IC phase. 
The exponential term represents a repulsion between do- 
main walls. Hence, only if the domain wall self-energy 
is negative the formation of domain walls is energetically 
favorable. The minimization of Eq. (6) with respect to n 
for a clamped crystal yields a continuous C-IC transition 
at 



a,iu) = -2ppl{u) = -2/3 



257(5 



2/(m-2) 



(11) 



in which the domain wall density vanishes as In ^ { [a{u) — 
ac{u)]/ac{u)}. 

Let us now turn to consider an undamped crystal. 
Its deformation is also a degree of freedom and, there- 
fore, we must minimize Eq. (6) with respect to u. In 
the C phase this deformation is Uc = ra/{2K(3), where 
K = K — r'^/{2j3). In the IC phase there is. in addition, a 
deformation e = u — Uc induced by the creation of domain 
walls. This is small close to the C-IC transition and so 
in Eq. (6) only lowest-order terms in e are relevant. 

Further results show that the strain dependence of the 
repulsion term can be neglected. Therefore, close to the 
C-IC transition it can be taken as 



AJ{u)n exp 



mp{u) 



2n 



~ a2nexp{—N/n), 



(12) 



where = rmrqo{uc) / 4: and 02 = —4Tra{uc)ac{uc)/ (mP). 
In contrast, the strain dependence of the domain- wall 
self-energy is essential close to the C-IC transition. This 
is such that 



(13) 



where Eq = £ [a{uc) - aduc)], Ei = £{r - igaduc)/ 
[{m - 2)a{uc)]} with £ = {2 - m)-Ka{uc) / {2m(3). As a 
result wc have that close to the C-IC transition Eq. (6) 
can be written as 

F ~Fe + y + {Eo + Ei€)n + a2nexp{-N/n), (14) 

where Fc = —a'^K/{4,K(i). Note that the strain depen- 
dence of the repulsion term yields an exponentially small 
correction to Ei. 

Minimizing Eq. (14) with respect to e one obtains 
that eeq = —nEi/K. In consequence, the free energy as 
a function of the domain- wall density has the form 

F =Fc + ain + a2ne^Tp{—N/n) — a^n^, (15) 

where a\ = Eq, has been defined before and 03 = 
E1/{2K). The term —0371,^ describes an attractive inter- 
action between domain walls and leads to a discontinuity 
of the C-IC transition as we will show. 

Let us mention that within the constant-amplitude ap- 
proximation used, the equation Ku + r{p'^) = which 
follows from Eq. (3) in absence of external forces taking 



5 = is not satisfied. In order to have it satisfied, higher 
order corrections to in the so-called small anisotropy 
parameter (e^ in the following) must be taken into ac- 
count. However, these corrections lead to higher order 
terms in £m in the thermodynamic potential^^ and there- 
fore are not important for the attraction. 

The two types of coupling between the order param- 
eter and the strain considered in Eq. (3) contribute to 
this attractive interaction. The ratio between these two 
contributions in the coefficient 03 is 



iti.tt'^ g qo{uc) 



{m-2)r o-(uc) 2-''(m - 2) r 



(16) 



where the parameter Sm = —{rn-K'^ /2'^)[aQ/ac{uc)\ de- 
fines the condition of weak anisotropy in the order pa- 
rameter space. ^ Eq. (6) is valid by virtue of the smallness 
of this parameter. 

The typical periods of the IC structures are larger than 
the atomic distances dat- usually they are such that 
<lo{uc) ^ 10^^(i~(^. The ratio g/r can be roughly es- 
timated as dat from dimensional arguments. Therefore, 
the relative importance of the two elasticity- mediated in- 
teractions is determined by the ratio between two small 
parameters qo{uc)dat/£m- As a result, the striction- 
mediated interaction is the most important one if 10~^ < 
£m < 1- Recall that there are also systems without Lif- 
shitz invariant (type II) where the striction-mediated at- 
traction is unrivaled. 

For m = 3, 4, 6 we have 



£3 = 



45 
J 



7 



1 2 



£4 



1 



£6 



2(3 



1/2 



(17) 



The smallness of qo{uc) — (uc) / {46) , in atomic units, 
implies the smallness of a{uc) since S, which gives the 
curvature at the minimum of the soft branch at the N- 
IC transition, has usually normal values. Therefore, the 
condition of weak anisotropy is easily fulfilled for m = 6 
even if all other coeSicients have normal atomic values. 
However, it is not the case for m = 4,3 which need, in 
addition, a special smallness of 7. 

2. Real systems 

We shall show that the free energy of real systems in 
terms of the domain- wall density has the same form that 
Eq. (15) close to the C-IC transition. To minimize the 
thermodynamic potential Eq. (3) with respect to the 
elastic degrees of freedom we first divide Uij into spatially 
homogeneous and inhomogeneous parts. The minimiza- 
tion should be carried out separately for these two parts 
because they represent degrees of freedom of the system 
independent one of each other. The above-mentioned 
division can be presented as 



(18) 
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where Ui{k) are the Fourier components of the displace- 
ment vector. Because only variations along the z axis are 
allowed we have that k = (0, 0, fc^)- 

Let us consider first the case in which the striction is 
the only coupling between the order parameter and the 
strain, i.e., we take gtj = in Eq. (3c). As a result of 
the minimization of Eq. (3c) with respect to the strain 
one obtains 



F2 = --^rijTkiXi^Jo - 



(19) 



in the Fourier space, where jk^ represent the Fourier com- 
ponents of the function j = (? and A~^j is the ijkl com- 
ponent of the inverse of the tensor Ayfe; . In the real space 
this is 



' 33 

^3333 



' 33 



2A 



3333 



■(/)• (20) 



Note that the first term of this functional is formally iden- 
tical to Eq. (4) if, for that case of infinite shear modulus, 
we take g = 0. Note also that the contribution of the 
second term to the thermodynamic potential can be un- 
derstood as a renormalization of the coefficient (3 in Eq. 
(3b). As a consequence, the free energy of a real crystal 
in terms of their domain-wall density has the form Eq. 
(15) with the renormalized constants 

/3 ^ /3-ri3/(2A3333): (21) 

r^K ^ TijrkiX^l- 4^X3333- (22) 

If both striction and Lifshitz invariant mediated cou- 
plings are taken into account, in addition of these renor- 
malizations we have 



g'/K 

rg/K 



~ 533/ ''^3333, 
''33533/^3333- 



(23) 
(24) 



In this case, after the minimization over the strain one 
obtains also terms that renormalize the coefficients of the 
invariants p'^V z'f and p'^iV z^)^- These invariants were 
not considered previously in Eq. (3) because their role 
is not essential. Therefore these contributions can be ne- 
glected. 

Let us discuss an isotropic case in which the shear mod- 
ulus ^ of the system is finite. The elastic moduli of such a 
system are Xiju = [K - {2^i/'d)]5ij5ki+ iJ,{5ikSji + dadjk)- 
One can see that the above mentioned renormalization 
of the coefficients in Eq. (15) results in 



4 7' 



3K{K + 



3K {K + 



(25) 
(26) 
(27) 



If the striction-mcdiated contribution to the attraction 
term of Eq. (15) is neglected, i.e., one takes r = 0, the 
resulting attraction term coincides with the previously 
reported in Refs.^'^. 

3. Other contributions to the attractive interaction 

For IC ferroelectrics thermal fluctuations of domain 

walls induce macroscopic electric fields which generate a 
van der Waals-like attraction between them.** Its contri- 
bution to free energy of the system is 



Fvdw = ' 



3 T 



1 T 
167r 'p 



In 



{I < ho), 



{I > ho), 



(28) 



where I is distance between domain walls, hn — 
dat[To/{<d-T)]^/'^ and it has been taken m = 6 (K2Se04- 
type crystals). 

Let us estimate, at the C-IC transition, the ra- 
tio between this contribution and the elasticity-induced 
one considering the striction-mediated part only (the 
striction-mediated attraction and that due to the Lif- 
shitz invariant arc expected to be of the same order of 
magnitude for K2Se04 as we shall see below). 

Prom the last term of Eq. (15), the elasticity- induced 
attraction between domain walls at the C-IC transition 
point is 



"el 



-(27r) 



K 



(29) 



where AK/K ^ 10~^ 10~^ is the relative change 
of the bulk modulus between the normal and the in- 
commensurate phases, dat is the interatomic distance, 
Ao = 27r/(7o(Mc) ^ lO^fiat is the typical wavelength of the 
IC phases at the N-IC transition, A/oc ^ Aq is the typical 
wavelength of the IC phases at the C-IC one, and Tq is 
equal to Tat ~ 10^ lO^K for displacive systems and to 
Tioc for order-disorder ones. Thus, the above mentioned 
ratio is found to be 



(10^ 10) 



Tloc (Q - Tioc) 



(10-=- 10-^)?)2c 



{hoc < ho), 

{hoc > ho)- 

(30) 



Taking into account that typically O — Tioc — Tioc, one 
finds that the contribution of van der Waals-like attrac- 
tion would be greater than striction-mediated one for 
order-disorder systems (To ~ Tioc), but it turns out 
to be as much as one ordcr-of-magnitude smaller than 
the striction-mediated contribution for displacive sys- 
tems (To ~ T„t). 

We mention that there is as well a universal 
fluctuation-induced attraction between domain walls op- 
erative by virtue of the solid-state elasticity.^ However 
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its effects are orders of magnitude smaller than the 
elasticity-induced interaction considered in this paper 
both for order-disorder and displacive systems. 

III. THERMAL ANOMALIES CLOSE TO THE 
C-IC TRANSITION 

A. Theoretical formulaes 

In this section we deduce formulas for thermal anoma- 
lies that a free energy like Eq. (15) describes. The pro- 
cedure is similar to that exposed in Ref.^^ but keeping in 
mind the different form of the attraction term. 

The coefficient ai of Eq. (15) can be expressed as 
fli = aiT{T - T,) where = 6 + [K / K)[a,{uc) / ar] 
is the temperature at which the domain- wall self-energy 
changes its sign. However, due to the domain-wall at- 
traction, Tc does not coincide with the C-IC transition 
temperature that will be determined below. 

The minima (stable or metastable) of the free energy 
Eq. (15) are obtained at certain values of the domain 
wall density, say, rie, such that the conditions 



dF 
dn 



= ai+ a2(l + Ar/ne)e-^/"= - 2a3ne = 0, (31) 
_ 2a3ne {N'^ - Nn^ - nf) - mN'^ 



are fulfilled. 



nl {N + Tie) 



Note that at Ti 



it 



IS 



>0, 
(32) 

Tc - 
obtained 



{d'^F/dr 



= and, therefore, it represents the limit 



of undercooling of the IC phase: below this temperature 
any minimum of the free energy (stable or metastable) 
is an state without domain walls. 

From Eq. (31) and the condition of continuity of the 
free energy at the C-IC transition point, i.e., F = F^, we 
obtain that 

Tioc = Tc - aznioc{N - nioc) / {air N) , (33) 
nioc = Nln-^ [a2Ar/(a3nL)] , (34) 

give, respectively, the C-IC transition temperature and 
the density of domain walls at the C-IC transition point. 
As we see, as long as 03 is not zero the C-IC transition 
temperature does not coincide with Tc and at this tem- 
perature the domain-wall density is not zero. 
The latent heat of the C-IC phase transition is 



Q = TiocttiTnio 



(35) 



and the specific heat of the IC phase depends on the 
temperature as 



(36) 



with a Curie constant 

Ci=TaiT{ne/Nf{N + ne). 



(37) 



Note that Cp diverges at the limit of undercooling of the 
IC phase. We mention that within a frame of a continu- 
ous C-IC transition the anomaly of the specific heat has 
the same form.^'^ 

For IC ferroelectrics, e.g., K2Se04-type crystals, the 
electric susceptibility in the IC phase is also interesting. 
The attraction term in Eq. (15) is proportional to the 
square of the domain walls concentration. Therefore, it 
does not contribute to the electric susceptibility since the 
electric field changes the distance between domain walls 
but not its concentration. Thus, from previous works (see 
e.g., Refs."'^^) and from Eq. (31) wc have approximately 
(taking He — nioc) the electric susceptibility as 



X 



C2 



T-T2 



where 



C2 = ^{N + nioc), 



To=Tc- 



aiT 



<Ti, 



(38) 

(39) 
(40) 



with Ps being the polarization within the domains. 

Note that at T = Ti, the temperature where the spe- 
cific heat formally diverges, the electric susceptibility re- 
mains finite. This might explain the experimental obser- 
vation reported in Rcf.^^, according to which the maxi- 
mum of specific heat occurs at a higher temperature than 
the maximum of the dielectric constant. 

It should be mentioned that the compressibility and 
the thermal-expansion coefficient also exhibit an anomaly 
similar to that of the specific heat. In other words, the 
compressibility of the crystal in the IC phase close to 
the lock-in transition can be quite high. However this 
is valid only for the static compressibility because the 
"softness" is associated with changes in the domain-wall 
lattice under the influence of the external pressure. Since 
the domain-wall processes are very slow one has to use 
very low frequencies measurements to probe this "soft- 
ness." 



B. Comparison with experiments: K2Se04 and 
Rb2ZnCl4 

There is an extensive literature devoted to the experi- 
mental study of the anomalies associated with C-IC tran- 
sitions (see, e.g., Ref.^^). Those which are focused on 
K2Se04 result of special interest for us due to the fol- 
lowing reasons. First, it is known that the free-energy 
expansion results in the form of Eq. (3) (with only one 
temperature-dependent coefficient) is valid to describe 
both the C-IC transition and the N-IC one.^"^ It allows 
us, in principle, to calculate the coefficients of the free 
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energy Eq. (15) independently on data of the C-IC tran- 
sition, i.e., the coeiEcients can be determined from the 
experimental data about the N-IC transition. Second, so 
far reported data permit to do it. In fact this valuable 
work was done previously by Sannikov and Golovko.^^ 
All coefficients of the functional Eq. (3) were consistently 
determined with the exception of gij , which needs more 
complete studies on the pressure dependence of the N-IC 
phase transition. According to the existing studies^** one 
can see that the above estimation g/r ~ dat is correct 
by an order of magnitude. Using the value sq ~ 10~^ 
given in Ref.^^ one sees from Eq. (16) that contributions 
of the two elasticity-induced mechanisms are compara- 
ble for K2Se04. Still in what follows we shall consider 
the striction-mediated interaction as the only one mecha- 
nism responsible of the attraction between domain walls, 
discussing at the end possible consequences of this as- 
sumption. Thus, our comparison between theory and 
experiments for K2Se04 can be made with no fitted pa- 
rameters and we know no other example where this is 
possible. 

It has been already mentioned that in the interpreta- 
tion of the experimental data an important characteristic 
of the "domain- wall systems" must be taken into account: 
they relax to equilibrium so slowly that it is impossible 
to observe their equilibrium states. It is evidenced in 
the hysteresis phenomena typical of the IC phases. It 
impUes that the observed thermal anomalies really cor- 
respond to smeared phase transitions predicted theoret- 
ically for systems in their equilibrium states. Hence, the 
aim of any comparison between theory and experiments 
must be to see if the experimental data can be reason- 
ably well explained as a result of the smearing of the 
equilibrium (theoretical) phase transitions. This fact re- 
solves the apparent contradiction between reported ex- 
periments so far, which show a great number of contin- 
uous C-IC anomalies. and the theory, which predicts 
discontinuous C-IC transitions. It motivates also a pref- 
erential comparison between integral characteristics such 
as the enthalpy of the transition. 

For K2Se04 whose non-zero elastic moduli are^^ 

Xii = {5.3,5.0,3.6,0.48,1.5,1.6}, 

(41) 

Al2,13,23 = {1.7,1.5,2.0}, 

in units of 10^^ dyn . cm~^, 

n/ao = {0.7, 2, -8, 0, 0, 0} • 10^ (42) 

(abbreviated notation is used and the index i runs the 
values 1,2,... 6) and (5K/{Kao) = 1.3- 10"^ dyn"! -cm^, 
we find that a3N/a2 ~ 10"^ xhen Eq. (34) gives 
n-ioc ~ O.IA^ and therefore the change of the "wave naum- 
ber" of the IC modulation would be such that qo /qioc ~ 2 
{qioc = 27rn;o(,/™)- From reported data in Ref.^''', the 
experimental value of this ratio is qq j qioc ~ 

3.5. The 

difference between both values can be commented in the 



following way. First note that the separation between 
the nearest points to the C-IC discontinuity in the ex- 
perimental curve of q{T) (Fig. 4 in Ref.^^) is about 
4 K. But as it can be seen in Ref.^*, all the observ- 
able C-IC thermal anomaly occurs into a region of about 
1 K. Therefore, one cannot be sure that the experimen- 
tal value really corresponds to the wave vector at the 
lock-in transition. It might be even more important the 
fact that the observed wave vector of the IC modula- 
tion never corresponds to its equilibrium value close to 
the C-IC transition because the relaxation times of the 
system are much greater than the measurement times. 
Thus, what is really meant by the experimental ratio 
Qo/Qioc remains somewhat undetermined. Second, it is 
quite possible that striction and Lifshitz invariant medi- 
ated contributions to as compensate each other. Indeed 
both transition temperatures and qo diminish with the 
pressure,^** a behavior that for an isotropic case evidence 
opposite signs of r and g. 

To the best of our knowledge, the most detailed study 
about the thermal anomalies of K2Se04 was made in 
Ref.^^. Unfortunately, the discussion of the reported data 
was based on a theory^^'^° in which the prediction of the 
discontinuity of the C-IC transition is erroneous because 
no discontinuity is obtainable if the considered model is 
treated consistently.^ We proceed to illustrate that the 
thermal anomaly for the specific heat of K2Se04 can be 
reasonably well understood as the smearing of a transi- 
tion which is discontinuous mainly due to the striction- 
mediated attraction between domain walls. 

For K2Se04 according to Ref.^'' it is obtained the value 
aiT — 0.05-^0.5 J-mol~^.K~^ and therefore, from Eq.(35) 
with Tioc = 95 K, the latent heat of the lock-in transi- 
tion is Q ^ 0.6 ^ 4 3 ■ mol"'^. This latent heat is enough 
to be observed experimentally if the transition (equilib- 
rium) were really observed, but its smearing is such that 
the thermal anomaly observed appears as continuous. 
In order to estimate the contribution of the anomaly of 
Cp to the enthalpy of the C-IC transition we calculate 
the difference between the C-IC transition temperature 
and the temperature at the limit of undercooling of the 
IC phase (the specific heat diverges at this limit). It is 
obtained a temperature difference of about 0.5 K. From 
Eq.(37) the Curie constant for the specific heat at the 
lock-in temperature is Ci ^ 0.1 ~ 0.03 J • mol~^. Bearing 
in mind that all the C-IC anomaly occurs into a region 
of about 1 K, the contribution of the Curie anomaly to 
the enthalpy results one order-of-magnitude lower than 
that of the latent heat. Therefore, theoretical enthalpy 
of the C-IC transition is in reasonable agreement with 
their experimental value^* 1.1 J • mol~^, and it can be 
said that almost all contribution to the thermal anomaly 
comes from the smearing of the latent heat. 

Another fairly- well studied compound is Rb2ZnCl4. 
However we do not dispose of the data necessary to cal- 
culate the coefficients of the free energy in the same way 
as the before for K2Se04. But we can use the experimen- 
tal value of the Curie constant C2 — 5 K, which seems 
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not to vary substantially for different samples (unlike to 
the Curie constant Ci, see e.g. Rcf.^^), as well as the 
data for = 0.12 iiC ■ cir-^ and nioc ^ 0.1 ^ 0.03iV (sec 
Refs.^^'^^) to estimate the values of the latent heat and 
the Curie constant Ci. Thus, from Eqs. (39), (37) and 
(35) with Tioc ~ 195 K, it is obtained Q = 7 ^ 2 J • mol"^ 
and Ci = 0.8 -j- 0.06 J • mol~^. These estimations are in 
a reasonable agreement with the experimental values re- 
ported in Refs.22,32 of Q = 6^2 J-mol'^ and Ci = 1^0.1 
J • mol~^. 

IV. CONCLUSIONS 

We have shown that the most universal interaction be- 
tween the order parameter and the elastic strain (stric- 
tion effect) leads to a long-range attraction between soli- 
tons. As a consequence, any C-IC transition in a solid 
should be discontinuous. This striction- mediated attrac- 
tion has been revealed in detail for a class of C-IC tran- 
sition that permits a clear and consistent analysis. The 
obtained results have been compared with experiments 



showing that the thermal anomalies in K2Se04 associ- 
ated to the C-IC transition can be well understood as a 
result of the smearing of this transition. 

We mention that we have considered the interaction 
between the domain walls in the bulk, neglecting the 
surface effects. The relaxation of the internal stresses 
and possible changes of the order-parameter distribution 
close to the surface can modify this attraction. The role 
of the striction effect in the surface properties of the IC 
phases deserves a special investigation. It seems promis- 
ing to study the role of the striction-mediated attraction 
in other systems, especially in the cases where contin- 
uous C-IC transitions are supposed for symmetries not 
admitting the Lifshitz invariant. 
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